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In many astrophysical problems, it is important to understand the behavior of functions that 
come from rotating (Kerr) black hole orbits. It can be particularly useful to work with the frequency 
domain representation of those functions, in order to bring out their harmonic dependence upon 
the fundamental orbital frequencies of Kerr black holes. Although, as has recently been shown by 
W. Schmidt, such a frequency domain representation must exist, the coupled nature of a black hole 
orbit's r and 6 motions makes it difficult to construct such a representation in practice. Combining 
Schmidt's description with a clever choice of timelike coordinate suggested by Y. Mino, we have 
developed a simple procedure that sidesteps this difficulty. One first Fourier expands all quantities 
using Mino's time parameter A. In particular, the observer's time t is decomposed with A. The 
frequency domain description is then built from the A-Fourier expansion and the expansion of t. We 
have found this procedure to be quite simple to implement, and to be applicable to a wide class of 
functionals. We test the procedure using a simple test function, and then apply it to a particularly 
interesting case, the Weyl curvature scalar xpi used in black hole perturbation theory. 

PACS numbers: 04.70.-s, 97.60. Lf 



I. INTRODUCTION 

The black holes which appear to exist in a wide range of masses throughout the universe (see, e.g. Refs. Q, H 
0, "l", 'H', "7, "s*!) are most hkely described by the Kerr solution of general relativity. The charged generalization is 
unlikely to be interesting, as macroscopic charged objects should be rapidly neutralized by astrophysical plasma. The 
Schwarzschild limit is an unrealistic idealization given how unlikely it is for an astrophysical macroscopic object to 
have precisely zero spin. This motivates a need to thoroughly understand phenomena in the vicinity of Kerr black 
holes. Such an understanding becomes quite important as studies probe ever more deeply into black holes' strong 
fields. 

Of particular interest to many applications is an understanding of Kerr black hole orbits. In the language of general 
relativity, "orbits" are bound, stable geodesic trajectories. It is a relatively simple matter to write down the equations 
governing these orbits and to integrate in the time domain to find the detailed trajectory that a body will follow. 

These orbits have a rich phenomenology, owing to the complicated shape of the hole's gravitational "potential" . At 
largish radii (r > 20 times the radius of the hole), a generic orbit is not too different from the ellipses of Newtonian 
theory. However, the plane in which this ellipse lies precesses (due largely to the spin of the black hole and the 
oblateness of the hole's geometry), and the ellipse precesses within that precessing plane. We can identify two 
fundamental orbital frequencies: a frequency 17^ characterizing the radial motion (from periapsis to apoapsis and 
back), and a frequency U,g characterizing the latitudinal motion. A third frequency of somewhat different nature 
describes the average secular accumulation of the angle about the hole's symmetry axis, and is denoted by O^. The 
various precessions of the orbit are due to mismatches between these frequencies: the orbital plane precesses at 
Vl^ — Vie; the orbital ellipse precesses at — $7^- Closed form expressions for all three of these frequencies have 
recently been worked out by W. Schmidt In the deep strong field of the hole, the frequencies become so different 
that the qualitative picture given above — a precessing ellipse on a precessing plane — ceases to be useful. The orbits 
just become complicated and messy. 

Despite this complicated nature, a wide class of functions of black hole orbits are completely described by the 
frequencies VLr and ^g. Any function of the form J[r{t),d{t)\ (a common functional form for black hole orbits, since 
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the metric is independent of both t and (jj) can be expanded as 

f[rit),9{t)] = J2 A„e-''="^'e-™"-* . (1.1) 

kn 

Unless otherwise noted, the index of all sums runs from —oo to oo. The fact that such expansions exist is very useful, 
since it suggests we can Fourier analyze a wide class of interesting orbit functionals to understand their harmonic 
dependence upon the orbital frequencies. 

Some functions have a more complicated form depending on all four components of the orbital worldline, = 
(t,r,9,(j)). A similar, but slightly modified, expansion can be written down which handles functions of this sort. 
Such an expansion is needed, for example, to give the harmonic decomposition of an orbiting body's stress-energy 
tensor, used in frequency domain perturbation theory of Kerr black holes |lOl |. One could also imagine using this 
harmonic expansion to describe the emission spectrum of hot material accreting onto a black hole. This could facilitate 
identifying features that are imprinted upon a black hole's x-ray spectrum. 

Actually computing the expansion coefhcients fkn turns out to be somewhat difficult. This is fundamentally because 
the r and 9 motions of a black hole orbit are coupled, and as a result are not periodic in coordinate time t (or proper 
time r). This difficulty can be fixed by working with a time variable A, recently suggested by Y. Mino jl3|, which 
decouples the r and 9 motions. With respect to A, the r and 9 motions are truly periodic. In contrast to the time 
domain expansion Hl.l|l . a similar expansion using A is straightforward to compute. 

The clocks of distant observers tick at evenly spaced intervals of t, not A. For the purpose of describing quantities 
that could be measured by such observers, the t expansion is more useful than the A expansion. Fortunately, it is 
straightforward to convert. That is the subject of this paper. The key observation is that observer time t contains 
oscillatory elements that are periodic with respect to Mino's time A. Thus, t itself can be expanded in a Fourier series 
of A-frequency harmonics. 

The remainder of this paper describes our prescription. In Sec. ^ we briefly discuss the i-domain description 
of the orbits. We then show how Mino's time A fixes many of the difficulties associated with these orbits in Sec. 
mil In Sec. IIVI we show how to use a A expansion to compute the the t expansion coefficients In Sec. we 
apply this technique first to a relatively simple function of black hole orbits, and then to the Weyl curvature scalar 
'04, demonstrating that everything works quite robustly. Appendix IXI discusses some important details related to 
implementation of these techniques. 

II. ORBITS IN BOYER-LINDQUIST TIME 

The geodesic equations that govern Kerr black hole orbits are usually presented in the following "classic" form '■ 
P'{^)' = [Eir^+a^)^aL,]' ^A[r^ + {L,-aE)' + Q]^Rir) , (2.1) 
p-ij^^^ ^ Q - cot^Ll ~ cos^ 9{1 - E'') = 0(9) , (2.2) 
m = csc^ 9L. + aE (^-^ - l) - ^ - '^(r^ ^) ^ (2-3) 



a sm 



-aL4l-'---^]^T{r,9). (2.4) 



Up to initial conditions, orbits are specified by the quantities E, L^, and Q ("energy", "z-component of angular 
momentum" , and "Carter constant" ) ; these quantities are conserved along any orbit of the family. For notational 
simplicity, we have put p"^ ~ + a? cos^ 9 and A = — 2Mr + a?. Note that Eqs. H2.1|) and H2.2|l have been divided 
by /i^, and Eqs. ^2.[\\ and (|2.4(l by p (where p, is the mass of a small body in an orbit); E, L^, and Q are thus the 
specific energy, angular momentum and Carter constant. The parameter r is proper time measured along the orbit; 
t is Boyer-Lindquist coordinate time. We choose < a < M; prograde and retrograde orbits are distinguished by an 
orbital inclination angle rather than the sign of the hole's spin. 

By picking initial conditions and physically reasonable values of the constants E, Lz, and Q, one can integrate these 
equations to obtain a worldline parameterized by proper time r along the orbit. Schmidt [9j has derived formulae 
for these constants as functions of an orbit's semi-latus rectum p, eccentricity e, and an inclination angle r, further 
discussion of these parameters is given in Appendix IXI Schmidt's formulae do not work well for circular orbits (e = 0). 
Formulae which apply to that case were originally worked out by Shakura ; we use a parameterization which was 
originally derived by Williams and then re-derived by Hughes [T^ . 



3 



For the purpose of understanding quantities which could be measured by distant observers, proper time is not 
a particularly good choice of parameterization for the orbit — it is connected to the orbit itself, and so contains 
components which oscillate with respect to the clocks of distant observers. Since the Boyer-Lindquist time coordinate 
t reduces at large radius to time as measured by distant observers, one should parameterize with t rather than t. It 
is trivial to convert: just divide the geodesic equations in t by dt/dr to obtain equations in t: 

dr dr f dt\ ^ 

di^d^[d^) ' ^^-^^ 

and likewise for dO/dt and d(j)/dt. Then, pick initial conditions and an allowed set of orbital constants {E, Lz, Q), and 
integrate to find z{t) = [r{t), 9{t), 

Using elegant Hamilton-Jacobi techniques, W. Schmidt |Qj has recently shown that bound orbits satisfying these 
equations are characterized by multiply-periodic motion in r, 0, and 4). These motions are given by three fundamental 
frequencies, fir, ^e, and 17^. In fact, the frequency Jl^ can be considered less fundamental than fij. and ^Ig. This 
is because the (j) orbital motion corresponds (in the language of Goldstein to a rotation-tyjie periodic motion, 
rather than an oscillatory or libration-tyjie periodicity. The frequency fi^ is the average rate at which accumulates 
over an orbit. Because d(j)/ dt depends only on r and 9, deviations from that average accumulation are oscillations at 
the r and 9 frequencies: 

4>{t) = n^t + Y, ¥'fc„e-*'="«*e-™"'-* . (2.6) 

kn 

Physically, one can imagine analyzing black hole orbits in a frame that co-rotates at the frequency fi^. In that 
corotating frame, the rotation-type periodicity at fl^ is removed, and only the libration-type oscillations at harmonics 
of fir and fig remain (see also discussion in Ref. pp. 466 - 467). 

By this logic, many functions f[z{t)] can be reduced to functions of r and 9 only. It is then possible to expand in 
a Fourier series as 

f[r{t), 9{t)] ^ /fc„e-^('="''+""'^)* . (2.7) 

kn 

Unfortunately, the functions r{t) and 9{t) are in general not periodic (although they are in the Newtonian limit where 
all the orbital frequencies are identical). This not-quite-periodic character is fundamentally due to the coupling the 
r and 9 motions in Eqs. H2.1|l and (|2.2|l : the functions {p^dt/dT)~^R and {p^dt/dT)~^Q each depend explicitly on 
both r and 9. (Note that this coupling remains if we use proper time along the orbit r as our parameterization.) The 
non-separated nature of the r and 9 motions makes it difficult to compute the coefficients f^n appearing in Eq. (|2.7|l . 
If the motions separated, one could define angl e variables w"^ = flrt and = flgt^ such that r would be a function 
only of and 9 a function only of ^2 ^2- Computing the coefficients fkn would then be straightforward (see, 
e.g., Ref. p. 466). Since the motions do not in fact separate, the angles w"^ and are not well defined. An 
alternative scheme to compute the Fourier series coefficients appears necessary. 



III. ORBITS IN MINO TIME 



In a recent paper, Y. Mino jl3j | introduced a new parameterization of Kerr geodesic motion which separates the 
and 9 motion. In terms of what we shall call "Mino time" A, the geodesic equations become 



dr^ ^ 



d9^ ^ 



i?(r). 


(3.1) 


9(0), 


(3.2) 


Hr,9) , 


(3.3) 


Tir,9) , 


(3.4) 



d\ 

dt 
d\ 

where R{r), 9(6'), <^ir,9), and T{r,9) are defined in Eqs. ||2IU-(E3I)- The r and 9 motions are now strictly periodic 
fimctions: 

r(A) = r(A-l-nA^) , 

9{X) = 9{X + nAg) , (3.5) 
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where n is any integer and the periods are given by 

The radial motion is taken to range between periapsis, Tpcri, and apoapsis, rap; the 6 motion ranges from a minimum 
^min to a maximum tt — ^min- (With a particular reparameterization, we can write the A^. integral in such a way that 
it behaves well as we approach the limit of circular orbits, rpori — > Tap. Likewise it is simple to reparameterize such 
that Ag is well behaved in the equatorial orbit limit, ^min tt/2. See Appendix IXI) 
For what follows, it will be useful to define the following frequencies conjugate to A: 

Tr,0 = 2^/A,,e, (3.8) 

as well as the angle variables 

w''''^ = Tr,eX ■ (3.9) 

These angles allow us to take advantage of the separated nature of r and 9 motion in Mino time: we treat r as a 
function only oi , 6 as a. function only of , and we treat w"^ and as independent parameters. This allows us 
to Fourier decompose any function of the orbital worldline using standard action-angle variable techniques |0 . 

Before moving on, we should analyze the remaining coordinate motions of black hole orbits — the observer (Boyer- 
Lindquist) time t and the azimuthal angle 0. Both of these motions consist of a component that accumulates secularly 
as a function of A, superposed on components which oscillate at Tj. and Tg. Let us analyze the oscillations first. From 
the geodesic equations and H3.4|l . we know that dt/dX and dip/dX are functions only of r and 9. This means that 
they can be expanded in a Fourier series: 

^ EE T(r, 9) Tfc„e-('=^^+"^'-)^ , (3.10) 



dX 

with the expansion coefficients given by 



^^Hr,9) = ^<i>,„e-('=^«+"^^)\ (3.11) 



kn 



Tkn = Tx^T^ rdw'- r dw'T[r{w'^),9{w')]e<''^'+^^'') , (3.12) 
(^"■j Jo Jo 

'fkn = -77^ r dw"- r dw' $ [r{wn,0{w')] e<'=«'''+"«'' ) . (3.13) 
(27r)^ Jo ^ 

In these equations and in what follows, r{w^) = r(A — /Tr) and 9{w^) = 9{X = /Tg). 
Because the functions T{r, 9) and <i>(r, 9) are real, we have the following relations: 

T_k,-n = fkn , (3.14) 

$_fe,_„ = $fc„ , (3.15) 

where the overbar denotes complex conjugation. The matrices Tfc„ and have another interesting property: Tko 
and Ton are non-zero, but Tkn = if /c 7^ and n 7^ (and likewise for <i>fe„). This lack of "crosstalk" between the 
9 and r harmonics is because T{r,9) and <i>(r, 0) have the form /(r) -|- g{9). To take advantage of this property, we 
define 

= ^feo , Tn = Ton ; (3.16) 
<fl EE $feo , K = *o„ . (3.17) 

Using the complex conjugate relations and Eqs. H3.16|l and (|3.17() . we rewrite the double sums appearing in the Fourier 
expansions (|3.11|) and (|3.1U|I as a pair of single sums jl8l |: 

,,00 00 
- ^ T{r, 9)=T + Y^ (T^-^^-^«^ + c.c.) -I- J] (T„^e--^-^ + c.c.) ; (3.18) 



dX 

k=l n=l 



00 



^ = $(r, 9) = T^ + J2 «e-'''^o^ + c.c.) + ^ (^^e"™^-^ + c.c.) . (3.19) 

fc=l n=l 
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The "c.c." means the complex conjugate of the preceding term. We have pulled the A; = 0, n = terms out of these 
sums and defined 

r - Too , (3.20) 

T0 = $00 . (3.21) 

These numbers tell us about the secular, average rate at which </> and t accumulate with respect to A. 
Using these results, it is simple to integrate for (f){\) and t{X): 

t{X) = TX + At{X) , (3.22) 

cl>{X) = T^A + A</)(A) . (3.23) 

We have chosen t{X = 0) = = <p{X = 0), and defined 

oo oo 

At{X) = ^ (At^e"'*^^^^ + c.c.) + ^ (Ai^e-'"^--^ + c.c.) ; (3.24) 

fc=l n=l 

oo oo 

A<^(A) = ^ {A4e-"'^o^ + c.c.) + J2 (ACe-'"''-^ + cc.) . (3.25) 

fc=l n=l 

We have defined Aff = iTj'" /{jTr,e) and A(j)'^f = i^y /{jTrfi). With this definition, we have separated t{X) and 
(^(A) into pieces which accumulate secularly with A plus pieces Af(A) and A(f){X) that oscillate at harmonics of Tg 
and Tr- 

Since f2<^ is the average rate at which (j) accumulates as a function of t and since F and are the average rates at 
which t and (f) accumulate as a functions of A, 

% = T^/r . (3.26) 

The other frequencies are likewise related: 

^0 = Te/r, (3.27) 

Vtr = T,./r . (3.28) 

When performing a harmonic decomposition of any function, we will want to work in terms of the angles = Tj A, 
for j = r, 0, 4>, and the average accumulated time T = FA. In terms of these variables, 

t{T,w\w') = T + At{w\w'), (3.29) 

(l){w'^,w'>,w'') = w't' + A(j){w^,w''), (3.30) 

where 

oo oc 

At{w', w^) = ^ (Atie-^'^^' + c.c.) + ^ (At^e-'"'"'^ + c.c.) ; , (3.31) 

k=l n^l 

DO OC 

A^{w',w^) = J2 (A^^e-^'^-' + c.c.) + J2 (ACe-"'"" + cc.) . (3.32) 

fe=l Tl=l 

Putting all of this together, the Fourier expansion coefficients fkn of any function of the form /[r(A), ^(A)] is 

fkn = -77^ r dw'' r d"^' M"^^)^ e{w')y^''^''+^^'') , (3.33) 

(27r)^ Jo Jo 

It is useful to note that the worldline z" can be reorganized in a similar form, separating the oscillations from the 
secular accumulations: 

z^{X) = zZ,{X) + A^«[r(A), ^(A)] , (3.34) 
where ^sec('^) = (r^^i 0, 0, T^A), and where 

Az"[r, 9] = [At{r, 0),r, 0, A0(r, 6)] (3.35) 
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can be expanded using the simple Fourier coefficients described by Eq. (|3.33|l with f = z". This leaves the worldline 
in the desirable form 

z"(A) = zZ,iX) + J2 AzLe'^('=^«+"^'-)^ . (3.36) 

kn 

By making use of Eq. H3.36|l . even rather complicated functional forms turn out to have a straightforward harmonic 
description. 

IV. CONVERTING FOURIER EXPANSION COEFFICIENTS 

There are two ways of Fourier expanding a function of the form f[r{t),9{t)] which are essentially equivalent: we 
can expand in observer time t, 

/[r(i),0W]=^/fc„e-*"-*; (4.1) 



or, we can expand in Mino time A, 



We have defined 



/[r(A),0(A)]=^A„e-^-^ (4.2) 



Tfe„ = kTg + flTr. (4.3) 

From the standpoint of measurable physics, the expansion (|4.1|l is more interesting — the components fkn tell us 
about the harmonic structure of / as seen by distant observers. However, the expansion H4.2|l is far more accessible 
— using Eq. (j3.33|l . it is straightforward to compute the expansion components fkn- In this section, we show how to 
convert the accessible components fkn into the measurable components fkn- 

We begin by taking the Fourier transform of f[r{t), 0{t)]. Using Eq. H4.1|l . we have 



1 f°° 

J2fknS{cj-iikn) = dtf[rit),eit)y^' 



1 r°° fit 

= ^ y ^rfA-/[r(A),0(A)]e-'(^). (4.4) 

Our goal is to evaluate the integral on the right-hand side of (14. 4|) and to find an expression relating fkn to fkn- To 
do so, we take advantage of the Fourier expansion for t{X) previously established, H3.22|l . 
We now insert dt/dX = T{r,9) and e^'^*'^) = e'"^^^ x e''^'^*(-^), under the integral: 

VA.„<5(c.-0fe„) - ^ rdA|T[r(A),0(A)]/[r(A),0(A)]e-^*(^)|e-rA^ 

1 r°° 

= — dXJ'[r{X),e{X),uj]e"^'^^ . (4.5) 

We next wish to insert the Fourier expansion of T[r{X),9{X),uj] under the integral. We first write this function in 
terms of the angle variables: 

T{w'',w',uj) = T[r{w'),0{w'')y^^*'''"'' '"""^ f[r{w''),0{w'')] . (4.6) 
The expansion of JF[r(A), 9{X),uj] is 

:F[riX),e{X),Lu] = J2 •^a6(^)e-*^-^ , (4.7) 
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where 

= TtAi r r dw' T{w\w\ uu) e--^^"-"- . (4.8) 
(271-)^ Jo Jq 

Inserting the expansion 1)4.6(1 into Eq. 14.5(1 . we find 

kn ab 

= :Fab{uj) d{ujr - Tab) , 

r-l^.Fafc(c^)(5(L^-17afc) , 

= r-i^.Fafc(r!„b)J(t^-r!„b) . (4.9) 

ab 

Equating the left-hand sides and right-hand sides, we read off 

/fcn-^fc„(r!fe„)/r. (4.10) 

We use this form of the Fourier expansion coefRcients in all of our calculations. 

V. EXAMPLES 

In this section we use the methods described above in two test cases. We first decompose and reconstruct a simple 
function of the form f{r,6). We then show how to decompose a more complicated function which appears in black 
hole perturbation calculations. 

A. Simple test case 

We now test our prescription by computing the expansion coefficients of a test function and showing that the 
reconstructed time series agrees with the original function. Our test function is C = cos 6, where r and 6 are 
the Boyer-Lindquist coordinates. Our calculations are actually performed using the tips given in Appendix In 
particular, we map the radial coordinate r to a coordinate ip defined via 

- (5.1) 



1 + £ COS -0 



a reparameterization commonly used in Newtonian orbital dynamics |20j . Whereas r oscillates from periapsis [rpori — 
pM/{l + e)] to apoapsis [rap = pM/{l — e)] and back, ip winds secularly from (periapsis) to tt (apoapsis) and beyond. 
We truncate all infinite sums at some finite value N, discussed below. We also map the 6 motion to a coordinate x 
via 



cos 



9 = \/ (a, E, l) cos x , (5.2) 



where z_ (a, E, i) is one root of a quadratic equation defined in Appendix^ The inclination angle l relates the angular 
momentum to the Carter constant Q: 

(5.3) 



Although L is not quite the geometrical angle describing an orbit's excursion from the equatorial plane, it is closely 
related, and has other convenient properties (see, e.g., Ref. for further discussion). 

It's worth noting that the parameterization 1(5.1(1 makes manifestly clear that the radial motion has a slowly 
converging Fourier expansion for large eccentricity. Using the binomial expansion. 



oo 



r = pM2^(-l)"£"cos" V • (5.4) 
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FIG. 1: The function ({t) — r{t) cos[9{t)] for orbits about a black hole with a = 0.9M, and with parameters p = 3, e = 0.2, 
L = 20°. Top panel shows this function computed directly from the geodesic equations (solid black line) as well as the time 
series constructed using Eq. 14.101 (dotted line; red in color version). In the reconstructed time series Crec(i), we have truncated 
the infinite sums ai N = 1. We find remarkably good agreement despite the small number of terms kept in the sum. The 
bottom panel shows the fractional residual, [("(t) — Croc(^)]/Cmax, for several values of A*'. The largest differences (red in color 
version) are for A'^ = 1, and have a magnitude of about 0.025. The next largest (green in color version) are for N — 2 and 
have a magnitude of about 0.003. The smallest differences (blue in color version) are for = 3 and have a magnitude of about 
0.0002. 



The amplitude of radial harmonic n is roughly e smaller than the amplitude of harmonic n — 1. When we truncate 
our sums at some finite value N, we expect that a reconstructed time series will have a fractional error of order e^''^'^ 
(the amplitude of the next neglected coefficient). Hence, many harmonics will be needed as e approaches 1. This slow 
convergence has been noted in studies of gravitational radiation reaction on eccentric black hole orbits . 

The top panel in Figs. Handle show the function ({t) computed in two different ways. The solid black line shows 
C(i) constructed by direct integration of the geodesic equations; the dotted line (red in the color version) shows ({t) 
reconstructed from a Fourier expansion using Eq. H4.1U|I . The lower panel of these figures shows the fractional residual, 
[C{t) - Croc(i)]/Cmax, whcrc Crcc(i) IS the reconstructed timeseries and (max = ?'max cos 6'min • 

Figure n compares C,{t) and Crcc(i) for an orbit with p = 3, e = 0.3, t = 20°; the black hole's spin parameter 
a = 0.9M. The reconstructed timeseries converges to (^(t) rather quickly: even when the infinite sums are truncated 
at iV = 1, the maximum deviation is only a few percent; the difference between ({t) and Cioc(i) for = 1 is barely 
discernible on the plot. The convergence is quite a bit slower when e = 0.6. The motion shown in Fig. |2]is for an 
orbit with p = 4, e = 0.6, and l — 50°. The timeseries is within 2% of ({t) when iV = 4; reducing the error by a 
further factor of ten requires increasing iV to 8. It's worth noting that, even when the amplitude error is relatively 
large, the phase alignment of CrGc(i) with ({t) appears to be very good. This bodes well for problems that rely on an 
accurate phase match between data and a model (or template) . 

Notice, in both Figs. ^ and ((2Jl, that the differences between ({t) and the reconstructed time series are noi precisely 
periodic: the wiggles in the lower panels of these figures do not repeat themselves regularly. This is a manifestation 
of the quasi-periodic nature of the orbital motion. As more terms are kept in these sums, we are more successful at 
capturing this quasi-periodic motion, and the magnitude of these wiggles quickly becomes small. 
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FIG. 2: The function ("(f) for orbits about a black hole with a = 0.9M and with parameters p = 4, e = 0.6, l = 50°. Top panel 
shows this function computed directly from the geodesic equations (solid black line) plus the time series constructed using Eq. 
14.1011 (dotted line; red in color version). We have truncated sums in the time series at A'^ = 4. We begin to see the need to 
keep a large number of terms at this relatively large eccentricity. The bottom panel shows the fractional residual for several 
values of N. The largest differences (red in color version) are for A*' = 4, and have a magnitude of about 0.02. The smaller 
differences (green in color version) are for N = 8 and have a magnitude of about 0.002. 

B. Black hole perturbations 

We now apply these techniques to a problem taken from black hole perturbation theory |0. Our goal is to 
understand how to decompose a complex function which describes how a small body perturbs the spacetime 
curvature of a Kerr black hole. From ip4, one can extract information about gravitational-wave emission and radiative 
backreaction on small compact objects orbiting massive black holes — extreme mass ratio binaries. When frequency 
domain perturbation theory is used to study this problem, tp^ is expanded in multipoles and in a harmonic series 
of the fundamental orbital frequencies. The gravitational waves generated by the system and their backreaction 
onto the orbit can then be extracted from that harmonic/multipolar expansion. Aside from being of great interest 
to the current authors, this example nicely illustrates the principles of this Fourier decomposition for functionals 
more complicated than the previous simple example. We will not dwell too much on the mathematics of black hole 
perturbation theory, but will point the reader to references where appropriate. 

The function ip4 can decomposed into multipoles as ^ .lOj 

Mtf,rf,ef,<l>f)=p-^Y1 J dojRi,n{rf,u:)Si,n{Of,aio)e'("'^f-'^'f'> , (5.5) 

Irn 

where the angular function 5* is a spin weighted spheroidal harmonic, and the radial function i? is a solution of a 
second order ordinary differential equation known as the Teukolsky equation (see Refs. [Toll2^ 12^ . I25ll2^ for 
a detailed discussion of the Teukolsky equation). The subscript / on the coordinates is a reminder that {tf,rf,9f,(j3f) 
denotes a field point. Coordinates without the subscript will refer to the location of a body orbiting the black hole. 
We will now rewrite ip4 completely in terms of sums, eliminating the need for the integral over to. 
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The radial function can be written in the form [22 1 123ll24 l25l l26| 

RUrf,u;) = Zliirf,^)R^irf,L,) + Zi^irf,u;)Rf^{rf,u;) , (5.6) 

where R^jl°° {rf,oj) are the two independent solutions to the source- free Teukolsky equation and where the functions 
Z are 

ZL{rf,io)^ fdt^^-'-^"^^''>hUr{t),0{t),rf,u;] (5.7) 



for ★ = H,oo. The function Ii^[r{t),9{t),rf,uj] depends upon the orbital worldline of the body perturbing the black 
hole spacetime. See Ref. |2^ for discussion in the case of a body in an equatorial, eccentric orbit; see Ref. 26] for the 
case of a body in an orbit that is inclined but of constant radius. (The general case, for orbits that are inclined and 
eccentric, is in preparation "27^.) 

We next rewrite Eq. ()5.7|l as an integral over A 



ZLirf,^) = j dAe^[-*(^)-™^(^)lir„Jr(A),0(A),r/,c.] , (5.8) 
where = I^^^dt/dX. Now we insert Eqs. (I3.22|) and 13.23|l into l|5.8|) so that we have 

ZLirf,^) = fdX e^^-^'-"^^*'hUrW,0{\),rf,Lu] , (5.9) 



where 

JUrW,dW,rf,io] =Xr„JKA),^(A),r/,He*^'^'^*[''(^)'^(^)l-"'^*['-(^)'''(^)l> , (5.10) 

with At and Acj) given by Eqs. 13.24|l and H3.25(l . Since J;*„ depends on A only through r(A) and 6'(A), it can be 
expanded as 

JL[r(A),0(A),r/,c.] =^ Jr^,„(r/,c.)e-^^-^ . (5.11) 

kn 

Putting this into Eq. H5.9|l and performing the integral gives 

ZLi^f,^) = 27r^ J/;„fe„(r/,w)5(u;r- T,„fe„) 

= "^i™fc«('^/,^)'5(w " n^nkn) , (5.12) 

jk 

where 

i^mkn = '^mkn/^ = + kfle + nflr ■ (5.13) 

Finally, when we substitute Eq. H5.12|l into Eq. (|5.5|l we obtain 

Mtf, rf,0f, cl>f) ^^Y. Ri^kn{rf)Si^kn{ef)e^^^'-''"^'"^'i'^ , (5.14) 



P 

Imkn 



where 



Slmkn{Sf) — Slm{6 f , a^mkn) 

Rhnkn{rf) = Z^^k^{rf)R'^{rf,n,nkn) + Z^t,^{rf)R^^{rf,VL„,kn) 

)/r . (5.15) 
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VI. CONCLUSION 



With the techniques described in this paper, it should now be a relatively simple matter to describe functions of Kerr 
black hole orbits in the frequency domain. Although the motion is not truly periodic with respect to observer time t, 
it is periodic with respect to Mino time A. It is thus quite simple to represent functions using frequencies T conjugate 
to Mino time. By using the fact that observer time t is itself periodic with respect to Mino time (after subtracting the 
secularly growing contribution), it is straightforward to convert the A-Fourier expansion into a t- Fourier expansion. 

As discussed in the Introduction, these techniques could find useful application to a variety of astrophysical problems 
involving Kerr black holes. One that is of particular interest to us is the problem of describing gravitational-wave 
emission from extreme mass ratio binaries p3 . 123 . I25I |2^ . Such systems are expected to be observable for 
future space based gravitational wave detectors. It should now be fairly straightforward to extend current black 
hole perturbation theory codes to handle the very interesting case of generic orbits — binaries in which the small 
body has both inclination with respect to the equatorial plane and non-zero eccentricity . If we had been unable 
to exploit the discrete harmonic structure of these systems, such a generalization would have had an enormous 
computational cost. Combining that analysis with a scheme to compute the evolution of the Carter constant (using 
a rigorous computation of a self force perhaps using a cruder approximation H^), it should then be possible 

to construct, in the adiabatic limit, the inspiral worldlines and waveforms followed by bodies spiraling into massive 
black holes. 
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APPENDIX A: PRACTICAL EVALUATION OF THE KERR GEODESICS 



In this appendix, we present some useful tools for handling functions of Kerr geodesies. A difhculty often encountered 
in evaluating these orbital motions is due to the presence of turning points in the motion: as the radial motion 
approaches periapsis and apoapsis, cir/d"time" passes through zero and switches sign (regardless of which time variable 
one uses). As the derivative approaches zero, one typically finds in a numerical evaluation that small stepsizes are 
needed to resolve the changing derivative; precision can be badly degraded in this case. The 9 behavior exhibits 
similar behavior due to the turning points at ^niin a-nd tt — ^min- 

A simple way to solve this behavior is to work with a functional form that automatically builds in the correct 
behavior as the turning points are approached. We first describe the transformation used to describe the 9 motion. 
The core idea of this transformation has been known for quite some time (29^ . and has been used extensively in 
work on circular Kerr black hole orbits it turns out to be particularly simple to use when studying geodesies 
parameterized by Mino time. We then show a simple transformation that greatly simplifies the description of the 
radial motion. This transformation has also been used quite a bit in previous work ,22. ■23J , but is worth discussing 
in the context of the Mino-time parameterization. 



1. Motion in 6 



We begin transforming the 9 motion by first defining the variable z — cos^ 9. Equation H3.2|l becomes 



d9 _ ^ /z2[a^{l-E2)]-z[Q + ^ + a^{l-E^)]+i 



d\ \ 1 — z 



The plus sign corresponds to motion from ^niin to tt — ^min, and vice versa for the minus sign. We have defined 
P = 0^(1 — E'^); z± are the two roots of the quadratic in the top line of Eq. l|Aip . 
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We next define the variable x'- z = z- cos^ x- x varies from to 27r, 9 oscillates through its full range of motion, 
from ^?rnin to TT — ^min and back. Examining dz/dO and dz/dx we see that 



dx 

de 



1 - z 



1 - z 



o<x<^; 

TT < X < 27r , 



(A2) 



Combining Eqs. (|Aip and l|A2|l . we obtain the geodesic equation for 

(be 

dX 



V/3(^+ - z) 



= V/3(2+ - z- cos2 x) 

Using Eq. 1A3() . it is straightforward to find A for all x- First, define 

1 



Ao(x) 



- F{7t/2 - X, 



note that 



(A3) 



(A4) 



(A5) 



In these equations, the function F((p, k) is the incomplete elliptic integral of the first kind, and K{k) is the complete 
elliptic integral of the first kind (using the notation of ^] ) . Then, 

A(x) = Ao(x) 0<x<7r/2 
2 



2 



if(v/z-/z+)-Ao(7r-x) V2<X<^ 



4 



if(v/z-/z+) + Ao(x-7r) ^<x<37r/2 



if(v/z_/z+)-Ao(27r-x) 37r/2<x<27r; 



also 



Ae 



i^(yZ7^) 



(A6) 



(A7) 



This form of Kg is perfectly well behaved even for orbits that arc confined to the equatorial plane (^min — 7''/2); this 
is not the case for the original form (|3.7|l . 

By combining Eqs. (|3.8|) . HA3|I . and HA7|I it is trivial to change variables so that integrals of become integrals 
over X'- 



w\x) = TeA(x) ; 
dw'^ _ ^ dA 
dx ^ dx 



(A8) 



27r 



Ae -y/ /5(z+ — z_ cos^ x) 



2i^(y;r7i7) vi-(z_/z+)cos2x ■ 

Equations ljA8|l and ljA9|l are used in our applications to perform all integrals with respect to the angle variable w' 



(A9) 
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2. Motion in r 



We use a similar trick to simplify the radial motion. First, we reparameterize the instantaneous orbital radius as 

pM 



1 + £ COS V' 



(AlO) 



Such a reparameterization is commonly used to study Keplerian orbits in Newtonian theory ;20j ; though relativistic 
orbits are not closed ellipses, the form (jA10|) remains very useful. The parameter e can thus be interpreted as the 
eccentricity, if) as the orbital anomaly, and p as the semi-latus rectum. As ip varies from to tt, r varies from periapsis 
(closest approach) to apoapsis (furthest distance): 

(All) 



' ap 



l + e 
pM 

l-e 



(A12) 



To proceed, we must do some massaging of the function R{r) defined in Eq. (|2.1|) . It is a quartic function of r, and 
thus has 4 roots: 

R{r) = {E^ - + 2Mr^ + [a^{E'^ - I) - LI - Q]r^ + 2M[Q + [aE - L,f]r - a^Q 

= il-E^){n-r){r-r2){r-r3){r-n) . (A13) 

The second line of Eq. ljA13|l is written in a way that is manifestly positive for bound orbits {E < 1). The roots are 
ordered such that ri > r2 > > r^; bound motion occurs for ri > r > r2. From these definitions, it is clear that 
ri = Tap, and r2 = rpori- 

The radii and r4 do not correspond to turning points of the small body's motion, but of course still represent 
zeros of the function R. (In fact, is typically inside the event horizon; when Q = or a = 0, r4 = 0.) It turns out 
to be useful to remap these radii as follows: 

r, = f^^. (A14) 

This remapping is simply for mathematical convenience; the parameters andp4 have no particular physical meaning. 
It is now a simple matter to derive the geodesic equation for ■(/;: 

dip _ MVI - E^ [(p-ps) -£{p + P3 cosip)]^^"^ [{p-pi) + e{p - picosip)]^^"^ 

dA ~ 1 -£2 

= PW;- (A16) 
As with the x reparameterization of the 9 motion, it is straightforward to find A(i/') using Eq. (1X1611 : 

In our applications, we evaluate this integral numerically. It is possible that an analytic form could be found in terms 
of elliptic integrals (though it appears to require more algebraic fortitude than these authors could muster). In any 
practical application, it is unlikely that such a form will be more useful or accurate than a numerical evaluation of 
(TnTll . 
Note in particular that 



This form of A^ is well-behaved in the limit of circular orbits. 
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Finally, we use these results to convert integrals over w'^ into integrals over -0: combining Eqs. (|3.8|l . IIA16|) . and 
l|XTH|l . we have 

w'itP) = TrX{^) ; (A19) 
dw^ ^ dX 

dip ^ dip 

We use Eqs. (|A19(I and l|A20|l to perform all integrals with respect to w'^ . 
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